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ABSTRACT 

This  report  describes  the  technical  accomplishments 
achieved  during  a  one  year  research  program  on  the  magnetic  and 
magnetoelastic  properties  of  materials  containing  elements  of  the 
rare-earth  series.  A  primary  goal  of  this  program  was  the  develop¬ 
ment  of  highly  efficient  high-frequency  ultrasonic  transducers 
based  on  the  extremely  high  magnetostriction  found  in  most  rare- 
earth  materials.  Work  has  been  carried  out  on  the  propagation 
of  elastic  waves  in  holmium  and  terbium  single  crystals,  and  re¬ 
sults  have  been  obtained  for  the  velocity  and  attenuation  of 
such  waves  as  functions  of  the  temperature  and  the  strength  of 
an  applied  magnetic  field.  Single- crystal  transducers  of 
gadolinium  and  terbium  have  been  constructed,  and  elastic-wave 
generation  in  these  transducers  has  been  observed,  although  only 
preliminary  results  have  been  attained.  A  new  theory  of  magnetic 
and  magnetoelastic  effects  in  the  rare  earths  has  been  developed, 
and  it  promises  to  shed  new  light  on  these  phenomena.  A  facility 
for  the  growth  of  large  single  crystals  of  rare-earth  materials 
has  been  constructed,  making  possible  the  study  of  a  wide  range 


of  materials. 
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I .  INTRODUCTION 

This  report  describes  the  technical  achievements 
attained  in  the  course  of  a  research  program  entitled:  "Research 
on  Rare-Earth  Metals  and  Compounds  and  Development  of  Applications 
Based  on  Their  Magnetic  Properties."  This  program  was  sponsored 
by  the  Advanced  Research  projects  Agency,  Department  of  Defense 
(ARPA  Order  1685),  and  it  was  administered  by  the  Army  Missile 
Command  under  Contract  No.  DAAH01-71-C-0259.  The  contract  per¬ 
formance  period  was  2  November  1970  through  1  November  1971. 

This  program  is  being  continued  for  a  second  year  under  Contract 
No.  DAAH01-72-C-0286,  beginning  5  November,  1971. 

The  principal  emphasis  in  this  research  program  was 
directed  toward  the  exploitation  of  the  extremely  large  magneto¬ 
striction  of  many  of  the  rare-earth  elements  and  compounds  for 
the  generation  or  ultrasonic  elastic  waves  at  high  intensity 
and  very  high  frequency.  The  development  of  rare-earth  magneto- 
strictive  ultrasonic  transducers  capable  of  generating  extremely 
intense  high-frequency  elastic  waves  will  find  many  applications 
in  such  areas  as  the  nondestructive  testing  of  materials,  invest¬ 
igation  of  the  nonlinear  elastic  properties  of  solids,  microwave 
communications  and  signal-processing  systems,  and  acoustic  holo¬ 
graphy.  In  addition  to  this  primary  concern  of  the  research 
described  here,  the  general  studies  of  magnetoelast ic  phenomena 
in  the  rare  earths  carried  out  in  this  program  will  undoubtedly 
have  important  consequences  for  the  development  of  n ew  permanent- 
magnet  materials,  such  as  the  rare-earth  cobalt  compounds. 


2 . 

The  research  carried  out  during  the  first  six  months 
of  the  contract  period  ha-  been  previously  reported  in  the  Semi¬ 
annual  Technical  Report^-  submitted  1  June  1971  in  accordance 
with  contract  requirements.  Those  projects  which  had  been 
successfully  completed  during  the  period  covered  by  the  Semi- 
Annual  Report  will  be  mentioned  only  briefly  in  this  report,  in 
which  emphasis  is  given  to  projects  not  completed  during  the  first 
six  months  and  to  new  projects  initiated  during  the  second  half 
of  the  contract  period. 

The  research  described  in  this  report  may  be  divided 
into  four  major  areas,  and  the  accomplishments  achieved  during 
the  contract  period  may  be  briefly  summarized  as  follows: 

1.  The  Growth  of  Single  Crystals 

Apparatus  for  the  growth  of  single  crystals  of  the 
pure  rare-earth  elements,  rare-earth  alloys,  and  rare-earth- 
transition-group  intermetallic  compounds  was  designed,  constructed, 
and  successfully  operated.  This  apparatus,  which  utilizes  the 
induction-heated  floating-zone  technique  for  crystal  growth, 
was  successfully  employed  for  the  growth  of  large  crystals  of 
dysprosium.  The  quality  of  the  crystals  produced  with  this 
apparatus  is  comparable  to  that  of  the  best  crystal  previously 
available  to  this  program.  Dysprosium  was  chosen  as  the  prin¬ 
cipal  material  for  the  initial  effort  to  produce  single  crystals 
mainly  because  it  is  the  least  expensive  of  the  magnetically 
interesting  pure  rare  earths  (Gd,  .Ub,  Dy,  Ho,  Er,  Tm)  and  it 
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is  readily  available  in  the  form  of  rods  suitable  for  the  tech¬ 
nique  employed  here.  Dysprosium  exhibits  a  solid-state  phase 
transition  (bcc-hcp)  at  a  temperature  slightly  below  its  melting 
point  and  has  a  rather  high  vapor  pressure  at  the  melting  point, 
two  properties  which  lead  to  most  of  the  difficulties  encountered 
in  attempts  to  grow  single  crystals  of  the  rare  earths.  The 
success  which  has  been  obtained  in  the  case  of  dysprosium  is 
extremely  encouraging,  and  it  is  believed  that  the  technique 
which  has  been  developed  in  this  program  can  now  be  successfully 
applied  to  the  remaining  magnetic  rare  earths  and  to  alloys  of 
these  elements  with  each  other.  No  attempt  has  yet  neen  made 
to  grow  single  crystals  of  the  more  complicated  intermetallic 
compounds,  although  preparations  have  been  started  to  attempt 
the  growth  of  crystals  of  the  extremely  interesting  compound 
TbFe2. 


2 .  Elastic-Wave  propagation  in  the  Rare  Sarths 

The  studies  of  elastic-wave  propagation  in  the  pure 
rare  earths,  which  were  reported  in  detail  for  the  case  of  terbium 
in  the  Semi-Annual  Report1  ,  have  been  continued  and  extended  to 
include  holmium  single  crystals.  The  measurement  of  elastic- 
wave  attenuation  and  velocity  in  the  range  of  several  MHz  is 
very  important  to  the  primary  goal  of  the  program,  since  the 
study  of  the  effects  of  applied  magnetic  fields  on  these  quanti¬ 
ties  leads  to  a  better  understanding  of  dynamic  magnetoelastic 
interactions  in  the  rare  earths.  Furthermore,  the  results 
presented  here  represent  the  first  comprehensive  study  of  the 
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elastic  properties  of  these  materials, 

3 .  Macnetostrict ive  Generation  of  Ultrasonic 
Elastic  Waves 

Magnetostrictive  ultrasonic  transducers  consisting 
of  thin  single-crystal  disks  of  both  gadolinium  and  terbium 
have  been  successfully  used  for  the  generation  and  detection 
of  elastic  waves  in  the  MHz  frequency  range.  Because  single¬ 
crystal  specimens  of  the  pure  rare  earths  exhibit  much  greater 
magnetostriction  than  do  polycrystalline  samples,  it  is  expected 
that  they  will  ultimately  perform  much  better  as  transducers 
than  do  the  thin  polycrystalline  films  which  were  studied  earlier 
and  reported  upon  in  the  Semi-Annual  Report‘d.  Consequently,  no 
further  work  was  done  on  the  thin  film  transducers  during  the 
second  half  of  the  contract  period  so  that  the  maximum  effort 
could  be  directed  toward  the  problems  associated  with  the  devel¬ 
opment  of  single-crystal  transducers.  Although  the  generation 
of  elastic  waves  by  means  of  magnetostriction  was  quickly  demon¬ 
strated  in  the  single-crystal  gadolinium  and  terbium  samples,  a 
detailed  study  of  the  characteristics  of  the  generation  mechanism 
could  not  be  completed  during  the  contract  period.  The  major 
problem  preventing  the  acquisition  of  reliable  data  was  the  inad¬ 
equacy  of  the  ultrasonic  spectrometer  normally  used  to  generate 
elastic  waves  in  piezoelectric  transducers.  This  apparatus  has 
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required  major  modifications  to  permit  the  application  of  a 
suitable  driving  signal  to  the  very-low-impedance  circuit 
required  to  drive  magnetostrictive  transducers.  New  equipment 
ideally  suited  to  the  task  will,  however,  be  available  during  the 
second  year  of  this  program. 

4 .  Theory  of  Magnetic  and  Maanetoelastic  Effects 
in  the  Rare  Earths 

In  attempts  to  explain  observations  on  the  generation 
and  the  propagation  of  elastic  waves  in  the  pure  rare  earths, 
which  have  been  studied  in  various  phases  of  this  research  pro¬ 
gram,  it  became  obvious  that  much  of  the  theory  which  has  been 
applied  to  these  problems  in  the  past  is  inadequate  for  the 
explanation  of  most  of  the  observed  effects.  This  inadequacy 
stems  in  part  from  the  understandable  attempt  to  treat  crystal- 
field  anisotropy  and  magnetoelasticity  in  the  case  of  the  rare 
earths  by  methods  which  were  appropriate  only  to  the  iron-group 
magnetic  materials.  By  treating  the  crystal-field  contributions 
to  the  magnetoelastic  Hamiltonian  in  a  more  realistic  way  than 
the  usual  effective-field  approach,  it  was  possible  to  see  imme¬ 
diately  how  the  more  correct  approach  would  explain,  at  least 
qualitatively,  the  results  which  have  been  observed.  As  a  fur¬ 
ther  contribution  to  the  theory  of  magnetic  and  magnetoelastic 
effects  in  the  rare  earths,  a  new  means  of  treating  dynamic 
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problems,  such  as  magnetic  resonance  and  ultrasonic  generation, 
has  been  developed.  This  new  technique,  the  details  of  which 
have  net  been  completely  worked  out,  will  permit  the  interpre¬ 
tation  of  most  of  the  resui.ts  obtained  thus  far,  but,  more  impor 
tantly,  it  will  permit  a  much  more  intelligent  approach  to  the 
question  of  predicting  what  materials  and  what  methods  will 
yield  the  best  results  in  attempts  to  design  magnetic  or  magneto 
elastic  devices. 

In  the  following  sections  of  this  report,  the  four 
areas  outlined  above  will  be  discussed  fully,  and  the  results 
obtained  during  the  contract  period,  both  experimental  and 
theoretical,  will  be  described  in  detail.  The  final  section 
of  the  report  summarizes  the  results  of  the  entire  research 
program  during  the  one-year  period  of  the  contract,  and  it 
attempts  to  put  the  program  described  here  into  the  context  of 
the  overall  ARPA  Rare-Earth  Program  and  of  the  other  current 
research  on  the  magnetic  properties  of  the  rare  earths. 
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II.  THE  GROWTH  OF  SINGLE  CRYSTALS 

Although  the  growth  of  single  crystals  of  the  rare 
earths  and  their  alloys  and  compounds  is  not  directly  related 
to  the  primary  goals  of  this  research  program,  the  ability  to 
produce  a  supply  of  crystals  of  various  compositions  and  orienta¬ 
tions  is  essential  to  the  overall  achievement  of  those  goals 
within  a  reasonable  period  of  time.  Even  though  single  crystals 
of  the  rare  earths  and  some  of  their  alloys  are  commercially 
available,  delivery  times  can  be  quite  long,  and  the  cost  of  such 
crystals  is  often  very  high  compared  to  the  cost  of  high-purity 
polycrystalline  material.  Furthermore,  in  the  case  of  most  of 
the  interesting  intermetallic  compounds,  single  crystals  cannot 
be  obtained  at  all  commercially  at  the  present  time .  Since 
nearly  all  materials  containing  significant  amounts  of  rare- 
earth  elements  are  highly  anisotropic  in  their  magnetic  proper¬ 
ties,  a  fundamental  understanding  of  these  properties  can  usually 
be  obtained  only  through  studies  of  the  properties  of  single 
crystals.  Thus,  a  program  to  permit  the  growth  of  single  crystals 
of  rare-earth  materials  of  high  quality  was  undertaken  at  the 
beginning  of  this  research.  Initially,  only  the  pure  rare  earths 
were  employed  in  this  project,  since  they  are  readily  available 
in  high-purity  polycrystalline  form.  Success  has  been  achieved 
in  this  project  in  the  case  of  dysprosium,  and  work  has  begun  on 
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the  magnetostrictd vely  very  interesting  compound  TbFe2. 

The  growth  of  single  crystals  of  the  rare  earths  is 
difficult  for  several  principal  reasons: 

1.  The  rare  earths  are  all  highly  chemically  reactive. 
Most  crucible  materials  interact  strong.ly  with  the  rare  earths 

at  temperatures  in  the  vicinity  of  the  melting  point.  It  was 
decided  to  overcome  this  difficulty  by  the  use  of  the  zone¬ 
melting  technique  in  the  floacing-zone  configuration. 

2.  Several  of  the  rare  earths  exhibit  very  high  vapor 
pressure  at  temperatures  near  the  melting  point  (approximately 

1  Torr  in  the  case  of  dysprosium) .  This  high  vapor  pressure 
makes  growth  impossible  in  a  vacuum,  thus  ruling  out  the  electron- 
beam  technique.  In  the  present  case,  it  was  decided  to  employ 
rf  induction  melting. 

3.  A  most  serious  problem  is  the  presence  in  most  of 
the  rare  earths  of  allotropic  phase  transitions  at  temperatures 
slightly  below  the  melting  point.  As  is  usual  in  such  materials, 
the  existence  of  such  transitions  greatly  complicates  the  problem 
of  growing  single  crystals.  In  the  case  of  the  rare  earths, 
however,  the  proximity  of  the  transition  temperature  to  the  melting 
point,  together  with  the  large  temperature  gradient  which  can 

be  maintained  in  the  zone-melting  technique,  apparently  help  to 
overcome  the  aforementioned  difficulties. 
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As  mentioned  above,  it  was  decided  to  employ  a  floating- 

zone  rf-induction  melting  technique  for  the  growth  of  rare-earth 

crystals.  Strain-annealing  methods  have  been  employed  for  this 
2-4 

purpose  ,  but  it  is  usually  possible  to  obtain  only  crystals 
of  random  size  and  orientation  by  this  method.  Other  conventional 

5 

methods  of  crystal  growth  have  proved  unsuccessful  .  One  of  the 
advantages  of  the  zone-melting  method  is  that,  if  it  is  successr>‘l, 
single  crystals  of  arbitrary  axial  length  can  be  produced. 

Although  the  diameter  of  such  crystals  will  always  be  limited, 
the  shape  of  the  crystals  resulting  from  this  method  is  ideal 
for  many  types  of  magnetic  measurement.  Several  of  the  rare 
earths,  such  as  gadolinium  and  terbium,  exhibit  vapor  pressure 
at  the  melting  point  low  enough  to  permit  the  use  of  electron- 
beam  melting  in  a  high  vacuum  environment,  but  since  the  remaining 
rare  earths  all  have  rather  high  vapor  pressures,  it  was  decided 
to  employ  rf  induction  melting  in  a  purified  inert -gas  atmosphere. 

Previous  attempts  to  grow  single  crystals  of  the 
rare  earths  have  been  reported  by  Savitskii  et  al  .  For  the 
most  part,  they  employed  the  strain-anneal  method,  but  they 
were  also  successful  in  growing  single  crystals  of  yttrium  by 
the  zone-melting  technique.  Crystals  obtained  from  the  Metals 
Research,  Ltd.,  Cambridge,  England,  are  apparent xy  grown  by  a 
zone-melting  method,  but  no  details  are  available  concerning 
their  actual  technique  of  growth. 
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The  initial  attempt  at  Rowing  a  dysprosium  crystal 
in  this  research  program  yielded  a  rod  of  6-nvm  diameter  and 
50-ram  length  containing  several  sizable  single-crystal  grains, 
randomly  oriented  .  This  result  was  obtained  with  only  one 
pass  of  the  zone  thoughout  the  length  of  the  sample,  because  the 
instability  of  the  industrial  rf  generator  made  control  of  the 
zone  length  quite  difficul. «  It  was  impossible  in  this  initial 
attempt,  therefore,  to  determine  the  effects  of  multiple  passes 
or  of  the  rate  of  growth.  The  size  of  the  single-crystal  grains 
was,  however,  encouraging,  since  they  were  all  large  enough  to 
be  used  in  most  of  the  other  experiments  being  carried  out  in 
this  program.  It  was  decided,  therefore,  to  modify  the  rf 
generator  to  permit  feedback  control  of  the  temperature  of  the 
molten  zone,  in  order  that  multiple  passes  at  any  rate  could 
be  made. 

The  construction  of  the  mechanical  system  for  the 
control  of  the  zone  position  and  rotation  of  the  sample  during 
the  growth  process  was  rather  straightforward.  The  system  was 
built  within  a  vacuum  system  originally  used  for  thin-film 
evaporation.  Since  the  sample  rod  was  to  be  rotated  during 
growth  in  order  to  insure  uniformity  of  the  molten  zone,  it  was 
necessary  to  provide  for  motion  of  the  rf  induction  coil.  This 
coil  was  fabricated  in  the  form  of  a  three-turn  coil  wound  as 
a  flat  spiral,  with  the  first  turn  soldered  to  a  flat  disk.  The 
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disk,  through  which  the  sample  passed  during  zone  melting,  served 
to  concentrate  the  rf  magnetic  field  so  that  a  molten  zone  of 
minimum  length  could  be  produced.  The  induction  coil  could 
be  translated  along  the  sample  axis  at  speeds  up  to  100  cm /hr. 

The  coil-sample  arrangement  j  3  shown  in  Fig.  1.  The  apparatus 
was  enclosed  within  a  belljar  which  could  be  evacuated  to  high 
vacuum.  The  bed  1 jar  could  be  out gassed  by  means  of  a  heating 
mantle,  and,  after  overnight  pumping,  an  ultimate  vacuum  pressure 

_7 

of  1.0  x  10  could  normally  achieved. 

Because  of  the  high  vapor  pressure  of  most  of  the 
rare-earth  materials,  crystal  growth  was  carried  out  in  a  helium 
atmosphere.  After  overnight  bakeout  of  the  system  at  high  vacuum, 
a  thin  film  of  the  specimen  material  was  deposited  over  the 
belljar.  me  vacuum  pump  was  then  isolated  from  the  belljar, 
and  helium  gas  was  introduced  into  the  chamber.  This  gas  was 
purified  by  passing  it  over  a  heated  rare-earth  specimen  and  through 
a  liquid-nitrogen  trap.  The  final  pressure  in  the  belljar  prior 
to  commencement  of  the  crystal  growth  was  approximately  500  Torr. 

As  mentioned  above,  the  most  serious  problem  in  the 
actual  zone-melting  process  is  the  control  of  the  zone  temperature. 
This  problem  is  illustrated  in  Fig.  2,  in  which  a  dysprosium 
sample  which  wa*.  zone-meiced  with  poor  control  is  compared  with 
one  obtained  with  good  temperature  control.  Adequate  control 


of  the  rf  output  of  the  induction  generator  was  achieved  by 
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means  of  a  saturable  reactor  placed  in  series  with  the  primary 
side  of  the  high-voltage  transformer  of  the  gnerator.  A  small 
control  signal  permits  variation  of  the  output  power  from  zero 
to  its  maximum  value,  and  the  output  is  easily  stabilized  by 
means  of  a  feedback  signal  proportional  to  the  output  current. 

With  the  very  stable  output  power  obtained  in  this  manner  from 
the  gen  'rator,  a  uniform  zone  could  be  maintained  with  manual 
control  of  the  power  level. 

An  attempt  was  made  initially  to  fabricate  dysprosium 
rods  of  6-mm  diameter  and  10  to  15  cm  length  by  induction  melting 
of  high-purity  material  in  a  water-cooled  copper  crucible.  It 
proved,  however,  difficult  to  obtain  reds  of  uniform  cross  section 
by  this  means,  and  this  method  was,  therefore,  abandoned.  Uniform 
cast  rods  of  vacuum-distilled  dysprosium  (and  of  the  other  mag¬ 
netic  rare  earths)  were  obtained  from  Research  Chemicals,  Inc., 
Phoenix,  Arizona.  These  rods,  6  mm  in  diameter  and  10  cm  in 
length,  were  cast  in  a  tantalum  crucible,  and  the  impurity  levels 
quoted  by  the  vendor  were  the  following: 


Y 

0.01  % 

Mg 

0.005  % 

Gd 

0.01  % 

Ta  0.3  % 

Tm 

0.01  % 

Si 

0.001  % 

F.r 

0.01  % 

Gaseous  Ele¬ 

Yb 

0.001  % 

Fe 

0.005  % 

Ho 

0.01  % 

ments  Unknown 

Ca 

0.02  % 

Sm 

0.01  % 

Tb 

0.01  % 

Since  it  was  stated  by  the  vendor  that  most  of  the  tantalum  in 
the  rods  was  probably  on  the  surface,  the  rods  were  etched  and 
electropolished  in  a  solution  of  1»0  %  perchloric  acid  in  ethanol. 
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The  rods  were  then  used  in  the  zone-melting  apparatus, 
and,  because  of  their  uniform  cross  section  and  density,  it  was 
possible  to  obtain  fairly  uniform  melting  on  even  the  first  pass. 
Three  or.  four  passes  were  made  in  most  cases  at  a  speed  of  approx¬ 
imately  4  om/hr.  Upon  examination,  after  etching  the  samples 
in  the  perchloric  acid-methanol  solution,  several  grains  could 
be  seen  on  the  outer  surface,  typically  1-3  mm  in  width  and  3-15  mm 
in  length.  A  section  taken  from  the  center  of  the  zone-melted 
region,  revealed  typically,  however,  a  large  interior  single¬ 
crystal  region  (60%  of  the  cross-sectional  area)  with  small  grains 
around  the  circumference.  A  longitudinal  section  revealed  that 
the  large  central  single-crystal  region  extended  several  centi¬ 
meters  along  the  rod  axis. 

It  was  decided  to  anneal  the  specimen  containing  the 
largest  of  these  single-crystal  regions  in  a  sealed  tantalum 
crucible,  and  this  process  was  carried  out  at  a  temperature  of 
1300°  C  for  66  hr.  Subsequent  analysis  of  the  sample  revealed 
that  the  central  single-crystal  region  had  grown  to  occupy  the 
full  5-cm  length  of  the  rod,  although  there  were  still  several 
small  grains  around  the  circumference.  The  average  cross-sec¬ 
tional  fraction  occupied  by  the  single  crystal  was,  however , 
estimated  to  be  greater  than  75%  after  the  anneal. 

X-ray  analysis,  using  back-reflection  techniques, 
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showed  no  evidence  of  a  mosaic  structure  or  of  other  major  lattice 
defects.  The  overall  quality  of  the  X-ray  results  was  comparable 
to  that  of  the  best  previous  crystals  availabJe  in  this  laboratory. 

The  problems  associated  with  the  growth  of  single 
crystals  of  the  rare  earths  have  not  been  entirely  solved,  but 
the  success  which  has  been  achieved  during  the  contract  period 
is  encouraging.  Crystals  as  large  as  any  which  were  previously 
available  have  been  grown,  and  their  quality  is  more  than  adequate 
for  most  of  the  experiments  associated  with  this  research  program. 
Work  in  this  area  will  continue,  and  it  is  believed  that  a  reason¬ 
ably  complete  characterization  of  the  crystals  produced  in  this 
apparatus  will  be  obtained  shortly.  New  problems  will  undoubtedly 
be  encountered  in  the  case  of  the  intermetallic  compounds,  such 
as  TbFe2,  currently  the  most  interesting  of  these  compounds  from 
the  point  of  view  of  this  program,  but  work  is  in  progress  to 
attempt  the  growth  of  this  and  other  materials,  and  these  problims 
should  soon  be  revealed. 
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III.  ELASTIC-WAVE  PROPAGATION  IN  THE  RARE  EARTHS 

Studies  of  elastic-wave  propagation  in  rare-earth 
single  crystals  provide  considerable  insight  into  the  nature 
of  the  magnetoelastic  interaction  in  these  materials.  In  fact, 
the  study  of  the  velocity  and  attenuation  of  high-frequency 
elastic  waves,  and  particularly  the  study  of  the  dependence  of 
these  quantities  on  the  value  of  an  applit 1  magnetic  field, 
provides  information  having  a  direct  bearing  on  the  phenomena 
whose  study  is  one  of  the  primary  goals  of  this  program:  the 
magnetos trictive  generation  of  high-frequency  elastic  waves  '  . 
Furthermore,  the  knowledge  of  the  values  of  the  elastic  constants 
obtained  from  measurement  of  the  elastic-wave  velocities,  is 
important  not  only  in  the  analysis  of  magnetoelastic  properties 
of  the  rare  earths  but  also  for  other  technological  applications 
of  these  materials.  The  work  described  here  has  included  the 
determination  of  the  elastic  constants  of  the  elements  terbium 
and  holmium  as  well  as  the  determination  of  the  magnetic-field 
and  temperature  dependences  of  elastic-wave  attenuation  in  these 
elements . 

A.  The  Elastic  Constants  of  Terbium  and  Holmium 
1.  Background: 

The  heavy  rare-earth  metals  undergo  magnetic  phase 
transitions  at  temperatures  in  the  cryogenic  range,  with  the 


16. 


exception  of  gadolinium,  whose  Curie  temperature  is  293  K. 

In  terbium,  a  paramagnetic-helimagnetic  transition  occuxs  at 

229  K,  and  a  helimagnetic- ferromagnetic  transition  occurs  at 

9  10 

221  K.  The  neutron-diffraction  studies  of  Koehler  et  al  ' 
revealed  that  in  the  helimagnetic  state,  the  magnetization  in 
each  basal  plane  is  uniform,  and  that  the  magnetization  direction 
varies  from  plane  to  plane  in  a  helical  manner,  the  interlayer 
turn  angle  varying  from  20.5°  at  229  K  (the  Nefel  temperature, 

Tn)  to  18.5°  at  221  K  (the  Curie  temperature,  T  )  .  Below  T^,, 
terbium  is  a  simple  ferromagnet  with  the  magnetization  lying  in 
the  basal  plane.  These  magnetic  transitions,  particularly  the 
first-order  transition  at  T  ,  which  is  accompanied  by  a  magneto- 
strictive  structural  change,  can  be  expected  to  have  a  signifi¬ 
cant  effect  on  the  elastic  properties  of  the  material  ,  and, 
indeed,  anomalies  in  the  isotropic  elastic  moduli  of  both  terbium 
and  holmium  have  been  previously  reported  at  temperatures  near 
the  magnetic  phase  transitions^. 

Similar  to  the  case  of  terbium,  holmium  undergoes  two 
major  phase  transitions.  There  is  a  paramagnetic-helimagnetic 

transition  at  T  =  131  K,  with  a  helimagnetic-ferromagnetic 
N 

transition  at  20  K.  Again,  this  structure  was  established 

12  13 

by  neutron-diffraction  studies  '  ,  which  revealed  that  the 

ferromagnetic  state  has  a  conical  magnetization  structure,  with 
a  helical  interlayer  variation,  but  with  a  uniform  projection  of 
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of  the  magnetization  along  the  hexagonal  axis.  Prom  20  K  to 
35  K  the  helical  struc'.ure  is  highly  distorted,  and  the  appli¬ 
cation  of  a  strong  magnetic  field  in  the  helimagnetic  range 

leads  to  the  establishment  of  several  complex  ordering  structures. 

14 

Several  of  the  elastic  constants  of  dysprosium  and 
gadolinium^'  ^  have  been  reported,  but  a  complete  set  of  values 
for  the  elastic  constants  of  terbium  and  holmium  has  not  pre¬ 
viously  been  available.  Only  the  velocities  of  ultrasonic 

longitudinal  and  shear  waves  propagating  along  the  c-axis  hive 
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been  reported  in  connection  with  attenuation  measurements  ' 

2.  Sample  Preparation  and  Experimental  Procedure: 

Single  crystals  of  terbium  and  holmium  were  obtained 
from  Metals  Research,  Ltd.,  Cambridge,  England,  in  the  form  of 
cylinders  6  mm  in  diameter .  These  crystals  were  spark-cut  to 
a  length  of  approximately  6  mm  and  electro-polished  so  that  the 
faces  we re  parallel  within  a  fraction  of  one  degree.  Three 
samples  were  required  of  each  of  the  elements  for  this  investi¬ 
gation: 

(1)  a  crystal  having  two  parallel  surfaces  normal 
■o  the  c-axis  [0001] ; 

(2)  a  crystal  having  two  surfaces  normal  to  the 
b-axis  [101.:]; 

(3)  a  crystal  having  two  parallel  surfaces  at 

a  45°  angle  between  the  c-plane  and  the  b-plane. 
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The  samples  were  eleotropolished  to  remove  surface  layers,  and 
their  orientation  was  determined  by  means  of  Laue  back-reflec¬ 
tion  X-rays  to  within  2°  of  the  desired  orientation. 

Measurements  of  the  elastic-wave  velocities  were 

.  .  19 

carried  out  by  the  pulse-superposition  method  of  Me  Skirnin 
The  experimental  arrangement  used  in  this  method  is  shewn 
schematically  in  Fig.  3.  Pulses  of  approximately  l~|jsec  dura¬ 
tion  and  of  varying  pulse  repetition  rate  were  generated  with 
a  tone-burst  generator  and  applied  to  the  quartz  transducer. 

The  transducers  used  in  this  work  included  both  X-cur  and  Y-cut 
quartz  crystals  with  a  20-MHz  fundamental  resonance.  These 
transducers  were  supplied  by  the  Valpey  Corp.,  Holliston,  Mass. 

The  transducers  were  bonded  to  the  experimental  samples  with 
Dow-Corning  276-V9  fluid.  In  cases  '•’•hen  the  attenuation  was 
quite  large,  other  bonds  were  tried,  but  they  failed  to  yield 
improved  results.  The  reflected  echoes  were  amplified  and  displayed 
without  rectification  on  an  oscilloscope,  and  the  repetition  rate 
was  adjusted  until  the  echoes  from  successive  pulses  were  in 
phase.  This  technique  permits  extremejy  accurate  measurement 
of  the  wave  velocity. 

rne  relationships  between  the  elastic  constants  and 

the  wave  velocities  for  all  the  crystallographic  orientations 

employed  here  are  given  in  Table  I.  For  the  calculation  of  these 

.  .  20; 

constants,  the  following  densities  were  used 


19. 


3  3 

Tbs  p  =  8.272  g/cm  ,  and  Ho :  p  =  8.803  g/cm  .  These  densities 
were  not  corrected  for  temperature  variations  in  view  of  the 
relatively  low  thermal -expansion  coefficients  of  these  metals. 

3.  Results  and  Discussion 

The  elastic  constants  obtained  in  this  investigation 

are  compiled  in  Table  II  for  terbium  and  in  Table  III  for 

holmium.  The  constants  c^3  and  c^  were  determined  directly 

from  measurements  in  the  [0001]  crystals.  The  constant  c  ^  was 

obtained  directly  from  measurements  in  the  [1010]  crystal,  but  in 

the  case  of  terbium  the  attenuation  of  the  longitudinal  waves  in 

this  crystal  was  so  high  at  temperatures  below  T  that  sufficient- 

C 

ly  strong  echoes  were  not  observable.  The  values  for  c^  in 

terbium  between  T  and  150  K  were,  therefore,  calculated  from  the 

measured  quasilongitudinal  and  quasishear  velocities,  v  and  v 

(cf .  Table  I)  ,  using  also  the  directly  determined  values  of  c^ 

and  c  .  In  both  holmium  and  terbium,  the  velocity  v  ,  together 

with  the  value  determined  for  c.^,  permitted  the  determination 

of  The  proper  combination  of  v^L  and  vqS>  together  v/ith 

utilization  of  c  „  and  cAA  yielded  the  value  of  the  elastic  constant 

33  44  ' 

c^  .  Tables  II  and  III  also  tabulate  the  normal  compressibilities 
of  the  two  elements.  The  formula  used  in  the  calculation  of 
these  compressibilities  is  given  in  a  footnote  to  Table  II. 

The  tabulated  results  for  terbium  are  also  shown  in 
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TABLE  I .  Relation  between  Measured  Sound  Velocities  and 
Elastic  Constants. 


Velocity 

Crystal  Axis 

Displacement 

Direction 

Relation 

(a) 

V1 

[0001] 

[0001] 

DVi2  .  c33 

(b) 

V2 

[0001] 

Any  direction 
in  the  (0001) 

=  C44 

(c) 

V3 

[  ioTo] 

[ioTo] 

pV32  =  C11 

(d) 

V4 

[  1010] 

[1210] 

*V4  '  C66  ’  ?<CU~  °12> 

(e) 

a 

> 

Quasilongitudinal  along  an 
axis  45°  between  [0001]  and 
[IOTO] 

pVQL  =  ^(C11+  °33  +  2  C44* 

+  4[i!cn'  c33>2+  (C13+C44,2lt 

(f) 

VQS 

Quasishear  along  an  axis  45° 
between  [0001]  and  [10101 

pVQS  “  *(cll+  C33  +  2  "44 ) 

-  ^i(=u-=33)2+(c13+=44)2)i 

TABLE  II.  The  Elastic  Constants  of  Terbium. 
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III.  The  Elastic  Constants  of  Holmium. 
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The  elastic  constants ^  (1)  through  (5),  are  given  in  units  of  101  dyne/cm  • 
compressibilities,  (6)  through  (8),  are  given  in  units  of  10'11  cm2/dyne. 
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graphic  form  in  Figs.  4  through  8.  In  the  paramagnetic  phase, 

the  elastic  constants  have  a  general  tendency  to  increase  with 

decreasing  temperature.  In  the  Ccise  of  c__  aid  c  .  a  sharp 

33  66 

minimum  is  observed  in  the  region  of  the  magnetic  transitions, 

14 

Figs.  4  and  5.  Rosen  and  Klimker  observed  a  similar  minimum 

at  T  for  c_„  in  dysprosium,  and  all  three  elastic  constants  in 
N  jj 

chromium  display  the  same  behavior  at  T  .  In  contrast,  the 

constants  c .  and  c  in  terbium,  shown  in  Fig.  7,  reach  a 
1^  13 

maximum  value  near  T^,  and  this  also  appears  to  be  the  case  for 
Cll'  Fig‘  6* 

The  error  in  the  constants  whose  values  were  obtained 
directly  from  the  measured  velocities  is  estimated  to  be  less 
than  1  %,  but  the  error  may  reach  3  %  in  the  cases  of  the  elastic 
constants  which  were  determined  indirectly.  The  special  dif¬ 
ficulty  in  the  determination  of  c^  directly  from  the  wave  velo¬ 
city  below  T  is  brought  out  in  Fig.  ,6.  The  difference  between 
the  results  for  c^  obtained  directly  and  indirectly  probably 
indicates  that  the  magnetoelastic  interaction  at  high  frequency 

depends  on  the  wave  propagation  and  polarization  vectors.  Such 

21 

dependence  was  observed  also  in  the  case  of  chromium  ,  where 
the  values  for  the  constant  c  determined  from  the  wave  velocity 
in  the  [100]  direction  differed  substantially  from  those  obtained 
by  measuring  the  velocity  in  the  [110]  direction. 


The  elastic-wave  velocities  in  terbium  measured  at 
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different  temperatures  are  given  in  Table  IV.  The  listed  values 
were  obtained  from  smooth  curves  of  the  velocities  as  functions 
of  the  temperature.  An  analysis  of  the  velocities  of  the  longi¬ 
tudinal  waves  along  the  c-axis  in  the  paramagnetic  phase  near 

T  reveals  that  they  follow  a  logarithmic  function  of  the 
N  ~ 

reduced  temperature,  as  shown  also  in  Fig.  9: 


fli 

v 

1 


-K  In  [ 


T—  T 


N 


N 


(1) 


This  rather  weak  temperature  dependence  is  in  agreement  with 

..  18 
theory 


The  observation  that  Eq.  (1)  is  obeyed  by  the  results 
is  taken  as  an  indication  for  the  predominance  of  a  single  relax¬ 
ation  process  in  the  vicinity  of  T^,  presumably  a  spin-phonon 
interaction.  The  relaxation  time  of  this  process,  T,  can  be 

computed  from  the  equation22 

v 

t  =  _L  .  ACC /(— 
ou2 

In  this  equation,  Av.^/v^  is  measured  at  T^,  and  Aa  denotes  the 
difference  in  the  attenuation  coefficient  between  the  paramag¬ 
netic  state  far  from  T  and  its  value  at  T  .  The  value  for 

N  S3 

A a  =  0.21  cm  \  determined  in  this  work,  leads  to  a  value  for 

S)  ^  i 

the  relaxation  rate  1/j  =  6  x  10  sec  ,  comparable  to  the 

9  —1  18 

rate  7  x  10  sec  obtained  by  Ltlthi  et  al  . 

The  temperature  dependences  of  the  elastic  constants 

of  holmium  are  shown  in  Figs.  10  through  13.  With  the  exception 
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TATI.K 

TV.  The 

HI trasonic 
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no 
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2.909 
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3.  58 
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1.692 
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150 
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1.675 
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1.735 

160 
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1.729 
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180 
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2 . 895 

1.716 

190 

2.989 

1.647 

1.662 

2.894 

1.712 

200 

2.968 

1.641 

1.654 

2.894 

1.708 

205 

2.958 

1 . 6  j8 

1.650 

2.894 

1.704 

210 

2.941 

1.635 

1.647 

2.894 

1.705 

216 

2.923 

1.632 

1.642 

2.893 

1.703 

218 

2.916 
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1.641 

2.893 

220 
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1.640 

2.893 

1.702 

221 

2.879 

222 
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1.632 

2.888 

1.640 

2.890 

224 

2.909 

1.630 

2.893 
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2.889 

226 

2.894 

1.629 

2.889 

1.635 

2.888 

1.700 

227 

2.882 

2.899 

1.625 

228 

2.895 

1.628 

2.892 

1.640 

2.887 

230 

2.912 

1.627 

2.896 

1.647 

2.887 

1.699 

231 

2.899 

235 

2.932 

1.625 

2.890 

1.648 

2.885 

1.701 

240 

2.942 

1.624 

2.891 

1.646 

2.884 

1.699 

245 

2.947 

1.623 

2.888 

1.645 

2.883 

1.698 

250 

2.950 

1.621 

2.884 

1.643 

2.881 

1.697 

260 

2.953 

1.618 

2.878 

1.639 

2.879 

1.693 

240 

2.954 

1.616 

2.874 

1.635 

2.876 

1,691 

280 

2.954 

1.613 

2.871 

1.631 

2.873 

1.688 

290 

2.955 

1.611 

2.866 

1.627 

2.870 

1.685 

300 

2.956 

1.609 

2.864 

1.623 

2.868 

1.683 

*  The  velocities  are  given  in  units  of  10^  cm/sec. 

For  explanation  of  the  symbols  used  for  the  velocities  see 
Table  I. 


°f  c13,  the  constants  all  exhibit  normal  behavior  in  the 
paramagnetic  phase.  As  is  approached,  the  constants  undergo 
changes,  the  largest  of  which  was  observed  in  c^,  Fig.  10.  In 
this  case  the  decrease  reaches  a  value  of  4  %  and  is  character¬ 
istic  of  a  second-order  phase  .transition. 

The  compressibilities  are  shown  in  Fig.  14  as  functions 
of  the  temperature.  These  compressibilities  were  calculated 
from  the  elastic-constant  values  given  in  Table  III.  It  is  note¬ 
worthy  that  the  changes  in  the  compressibilities  of  holmium  at 
T^  are  rather  small  compared  to  those  observed  in  terbium.  Fig  8. 
This  difference  can  be  attributed  to  a  large  variation  of  the 
compressibilities  near  the  Curie  temperature  and  a  relatively 

small  one  near  the  Nd'el  temperature.  This  behavior  was  actually 

14 

observed  in  the  case  of  dysprosium  by  Rosen  and  Klimker  .  The 
large  variation  in  holmium  at  T  was  not  observed  in  the  work 

w 

reported  here,  since  T  =  20K,  and  the  measurements  given  here 

v 

extends  only  as  low  as  78  K.  In  terbium,  however,  the  helimag- 
netic  phase  extends  over  a  temperature  range  of  only  8  degrees, 
so  that  the  effects  of  the  two  phase  transit iona  on  the  two 
compressibilities  overlap. 

A  knowledge  of  the  elastic  constants  makes  a  calcu¬ 


lation  of  the  Debye  temperature  feasible.  By  the  method  of 


direct  integration  for  hexagonal  crystals  ,  using  the  elastic 


constants  measured  at  300  K  listed  in  Table  II,  the  Debye 
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temperature  for  terbium  was  found  to  be  176.6  K.  This  is  in 

24 

good  agreement  with  the  value  174  K  computed  by  Rosen  from 

the  elastic  moduli  of  polycrystalld ne  terbium.  Both  these 

values  are,  however,  higher  than  the  calorimetrically  determined 
,  25 

value  of  158  K  .  Calculation  of  the  Debye  temperature  for 

holmiur.i  by  the  same  method  yielded  a  value  of  181  K,  which  is  in 

agreement  with  the  value  184  K  obtained  from  the  moduli  of  poly- 

24 

crystalline  holmium  .  As  in  the  case  of  terbium,  the  computed 

values  exceed  the  value  of  161  K  which  has  been  dtermined  from 

25 

specific-heat  data 

B.  The  Attenuation  of  Ultrasonic  Elastic  Waves 
in  Terbium  and  Holmium 

Changes  in  the  elastic-wave  attenuation  are  quite 
marked  near  the  magnetic  phase  transitions  in  the  heavy  rare 
earths.  The  variation  of  this  attenuation  with  the  intensity 
of  an  applied  magnetic  field  generates  valuable  information 
concerning  the  magnetoelastic  interactions  in  these  materials. 

In  fact,  the  attenuation  is  directly  related  to  those  magneto¬ 
elastic  properties  which  are  most  important  for  the  generation 
of  elastic  waves  through  the  magnetoelastic  interaction.  The 

results  of  tne  measurement  of  elastic-wave  attenuation  in  terbium 

1 

were  presented  in  detail  in  the  Semi-Annual  Report  .  Additional 
work  was  carried  out  in  this  area  during  the  second  half  of  the 
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contract  period,  and  measurements  on  the  attenuation  in  holmium 
were  begun.  These  results  are  summarized  here. 

1.  Effect  of  Magnetic  Field  on  Elastic-Wave 
Propagation  in  Terbium: 

In  order  to  measure  the  attenuation  of  elastic  waves 
in  terbium,  essentially  the  same  apparatus  was  used  as  that 
employed  in  the  measurement  of  velocity.  The  cryostat  was  in¬ 
stalled  in  such  a  way  that  the  sample  was  in  the  gap  of  a 
large  electromagnet  capable  of  producing  a  uniform  field  as 
high  as  17  kG.  The  measurements  on  terbium  were  made  in  the 
temperature  range  from  280  K  to  180  K,  which  includes  the 
paramagnetic,  helimagnetic,  and  ferromagnetic  phases  of  terbium. 
The  magnetic  field  was  always  parallel  to  the  basal  plane,  and 
it  was  parallel  to  either  the  a-axis  or  the  b-axis. 

Figures  15  and  16  show  the  temperature  dependence 
of  the  first  reflected  echo  amplitude  in  the  presence  of  a 
magnetic  field.  At  zero  magnetic  field,  the  echo  amplitude 
decreases  sharply  as  the  transition  into  the  helimagnetic  phase 
is  approached,  and  it  then  decreases  more  slowly  to  exhibit  a 
minimum  at  Tc,  with  a  shoulder  at  T  .  In  the  presence  of  a  mag¬ 
netic  field,  both  the  amplitude  of  the  minimum  and  its  position 
on  the  temperature  scale  are  affected.  At  small  values  of  the 
applied  field,  the  amplitude  of  the  minimum  increases  with 
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increasing  field  until  the  echoes  disappear  completely  at  a 
field  of  2  kOe  (the  field  is  always  given  in  terms  of  the 
applied  field,  since  the  sample  geometry  precluded  the  use  of 
a  simple  demagnetizing  factor  to  permit  the  calculation  of  the 
internal  field) .  An  increase  in  the  magnetic  field  beyond  this 
level  leads  to  a  decrease  in  the  amplitude  of  the  minimum.  The 
position  of  the  minimum  is  always  shifted  to  higher  temperature 
as  the  field  is  increased;  this  shift  is  largest  v/hen  the  field 
is  applied  along  the  a-axis. 

Figure  15  illustrates  the  dependence  of  the  temperature 
at  which  the  minimum  echo  amplitude  occurs  on  the  magnetic-field 
intensity  for  both  the  a-  and  b-axis  cases.  From  this  figure, 
it  can  be  seen  that  the  data  points  for  each  case  lie  on  a  straight 
line.  The  slope  of  the  line  representing  the  a-axis  data  is 
greater  than  that  for  the  b-axis  case.  The  intercept  of  both 
lines  with  the  temperature  scale  at  H  =  0  lies  at  T  =  228.5  K, 
approximately  the  value  of  T^.  The  difference  in  the  slope  of 
the  two  lines  is  undoubtedly  the  result  of  the  basal-plane 
magnetocrystalline  anisotropy . 

In  order  to  compare  the  results  shown  in  Figs,  15  and 
16  with  existing  theories,  the  attenuation,  rather  than  the  echo 
amplitude,  must  be  accurately  determined.  This  determination 
could  not  be  made  in  the  many  cases  where  only  one  echo  could  be 
observed  or,  of  course,  in  cases  where  the  attenuation  becomes 
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so  large  that  the  entire  echo  pattern  disappears .  In  the 
extreme  field  limits,  namely  H  =  0  and  H  =  15  kOe,  reproducible 
attenuation  measurements  could  be  made.  These  results  are  shown 
in  Fig.  18. 

Several  theories  have  been  developed  to  interpret 
elastic-wave  attenuation  at  magnetic  phase  transitions  .  All 

of  these  theories  are  based  on  a  critical  scattering  mechanism 
at  the  phase  transition,  and  it  should  be  observed  that  this  mech¬ 
anism  is  probably  inoperative  at  temperatures  well  removed  from 
the  phase  transitions  in  the  rare  earths ,  The  existing  theories 
all  lead  to  the  same  form  for  the  e  .t.ected  dependence  of  the 
attenuation  of  longitudinal  waves  on  wave  number,  q,  and  reduced 
temperature,  e  =  (T—T^J/T  ; 

a^q,  e)  =  Bq2e  ^  ,  (3) 

where  B  is  a  temperature-independent  constant  from  which  the 
magnetoelastic  coupling  can  be  determined,  and  ^  is  an  exponent 
which  characterizes  the  spin-fluctuation  mechanism  responsible 
for  the  magnetic  transition. 

In  order  to  compare  the  experimental  results  shown  in 
Fig.  18  with  Eq.  (3),  the  temperature-independent  contribution 
to  the  attenuation  was  considered  to  be  constant,  and  it  was 
subtracted  from  the  measured  values  to  give  the  contribution 
pertinent  to  the  theory.  Figure  18  gives  a  log-log  plot  of  the 
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attenuation  versus  T-T  ,  where  T  is  the  position  of  the 

P  P 

attenuation  peak  on  the  temperature  scale.  From  this  figure, 
it  can  be  seen  that,  in  the  absence  of  a  magnetic  field,  a  straight 
line  represents  a  good  fit  to  the  data.  In  the  case  of  an  appre¬ 
ciable  applied  field,  however,  two  power=law  regions  characterize 
the  results.  In  each  region,  the  exponent  is  independent  of  the 
direction  of  the  magnetic  field.  Furthermore,  the  exponent  of 
the  high-temperature  power  law  is  equal  to  the  exponent  which 
applies  at  zero  magnetic  field.  In  the  low-temperature  region, 
both  B  and  r>  ate  much  smaller  than  in  the  high-temperature  region. 
The  values  for  the  constants  B  and  rj  in  Equation  ( 3) 

obtained  in  this  work  are  compared  below  with  the  values  experi- 
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mentally  obtained  by  Pollina  and  LUthi  and  with  whose  calculated 
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by  Laramore  and  Kadanoff  : 

EXPONENT  n 

32  31 

Applied  Present  Investigation  P.  and  L.  '  L.  and  K. 

Field  High  Temp.  Low  Temp. 

0  1.415  ±  0.2  1.24  ±  0.1  3/2  or  4/3 

15  kOe  1.45  ±  0.2  0.32  ±  0.1 

CONSTANT  B 


0 

1.4  x  10  9 

1.2  x  10"9 

4.7  x  10  9 

15  kOe 
(a-axis) 

5.1  x  10“9 

2.3  x  10“7 

15  kOe 
(b-axis) 

10.0  x  10“9 

2.6  x  10"7 
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2.  Attenuation  of  Elastic  Waves  in  Holmium? 

Attenuation  measurements  in  holmium  have  not  yet  been 
completed.  The  attenuation  of  20  MHz  longitudinal  waves  propa¬ 
gating  along  the  hexagonal  axis  has  been  measured  in  the  temp¬ 
erature  range  from  300  K  to  78  K,  with  the  results  shown  in  Fig  20. 
It  can  be  seen  that  the  attenuation  decreases  slowly  with  decreasing 
temperature,  but  that  it  begins  to  rise  sharply  as  the  helimagnetic 
phase  transition  is  approached.  Below  T  ,  the  attenuation  again 

decreases  to  a  value  near  the  high-temperature  value.  This  peak 
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was  observed  previously  by  Pollina  and  Ltithi  at  frequencies 
from  50  MHz  to  170  MHz. 

The  effect  of  an  applied  magnetic  field  on  the  attenua¬ 
tion  of  c-axis  longitudinal  waves  in  holmium  was  also  investi¬ 
gated  in  the  temperature  range  from  4.2  K  to  300  K.  The  field 
was  applied  in  the  basal  plane  along  the  a-  or  b-axis.  In  the 
temperature  range  between  300  K  and  70  K,  where  quantitative 
measurements  we re  possible,  no  effect  of  the  magnetic  field  could 
be  observed  up  to  an  intensity  of  18  kOe,  nor  was  any  effect 
observed  down  to  35  K,  although  quantitative  measurements  were 
impossible  in  this  range.  This  behavior  is,  of  course,  contrary 
to  what  was  found  in  terbium,  where  the  effects  of  applied 
maqnetic  fields  are  very  strong,  even  at  temperatures  well  above 

the  transition  temperature,  T^« 

Below  35  K,  and,  in  particular,  below  20  K,  T  ,  the  effect 

C 
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of  the  applied  field  is  observable,  although  accurate  measure¬ 
ments  have  not  yet  been  possible.  Work  is  proceeding  in  this 
area,  and  it  is  hoped  that  some  light  can  be  shed  on  the  strik¬ 
ingly  different  behavior  of  holmium,  compared  to  terbium,  when 
measurements  are  made  with  different  crystal  orientations  and 
with  shear  as  well  as  longitudinal  waves. 


TRANSDUCERS 


In  the  Semi-Annual  Report  covering  the  first  six  months 
of  the  contract  period,  the  work  which  had  been  carried  out  on 
the  generation  of  elastic  waves  in  thin-film  polycrystalline 
magnetostrictive  transducers  .  of  various  rare  earths  was  reported 
in  detail o  The  emphasis  during  the  remainder  of  the  contract 
period  has  been  placed  on  the  generation  of  elastic  waves  in 
single-crystal  transducers 0  Although  the  single-crystal  trans¬ 
ducers  should  exhibit  higher  efficiencies  in  the  generation  of 
elastic  waves  at  high  frequencies  than  do  the  polycrystalline 
thin  films,  certain  experimental  problems  have  prevented  a 
thorough  study  of  this  generation  in  such  single-crystal  trans¬ 
ducers  during  the  period  of  interests  Elastic-wave  generation 
has  been  observed  at  frequencies  of  3-10  MHz  in  transducers  of 
gadolinium  and  terbium,  in  the  form  of  thin  single-crystal  disks, 
but  the  dependence  of  the  amplitude  of  the  generated  waves  on 
both  temperature  and  applied  magnetic  field  is  quite  complex,, 

Eoth  longitudinal  and  transverse  waves  are  generated  with  differ¬ 
ent  temperature  and  field  dependences*  Several  major  problems 
arose  during  the  measurement  of  the  elastic-wave  amplitudes  in 
this  work,  however,  and  the  data  obtained  aro  not  sufficient  at 
this  time  to  permit  a  satisfactory  interpretation*  The  results 
obtained  up  to  the  pres.mt  are,  therefore,  summarized  in  this 
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report/  together  with  a  description  of  the  improvements  which 
are  planned  to  the  experimental  technique ..  These  improvements 
will  make  both  the  acquisition  and  the  interpretation  of  data 
much  simpler o 

In  the  case  of  both  gadolinium  and  terbium,  elastic- 

wave  generation  has  been  observed  in  single-crystal  magnetostric- 

tive  transducers  at  temperatures  from  the  ordering  temperatures 

(T  =  293  K  for  Gd,  and  T^  =  229  K  for  Tb)  down  to  77  K,  at 

frequencies  in  the  range  of  severa?.  MHz„  Both  longitudinal  and 

transverse  waves  were  observed,  but  the  dependence  of  their 

amplitude  on  field  and  temperature  agrees  with  expectations  only 

in  the  case  of  terbium  in  the  ferromagnetic  phase „  In  this  case, 

the  generation  of  elastic  waves  exhibits  a  resonant  behavior 

with  a  resonant  value  of  the  magnetic  field  whose  temperature 

dependence  follows  the  same  type  of  curve  as  that,  observed  by 
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Blackstead  and  Donoho  for  ferromagnetic  resonance  xn  dyspro¬ 
sium.,  The  resonant  field  increases  with  decreasing  temperature 
below  the  Curie  temperature  until  it  reaches  a  maximum  value  at 
a  temperature  equal  to  approximately  two-thirds  the  Curie  temp¬ 
erature.,  The  value  of  the  resonant  field  then  remains  constant 
as  the  temperature  is  further  decreased »  This  behavior  is  obtained 
only  with  the  applied  field  along  a  hard  magnetization  direction 
in  the  basal  plane  (the  a-axis,  in  the  case  of  Tb)  and  with  the 
driving  rf  magnetic  field  along  an  easy  direction  (the  b-axis, 
in  the  case  of  Tb) .  This  behavior,  although  similar  to  the  case 


of  ferromagnetic  resonance  in  dysprosium,  is  totally  unpredicted 
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by  current  theory  ,  and  it  has  helped  to  stimulate  the  theoreti¬ 
cal  work  on  the  dynamic  magnetic  and  magnetoelastic  properties 
of  the  rare  earths  reported  in  the  following  section  of  this 
report . 

The  work  in  this  area  has  been  hampered  by  the  nature 
of  the  apparatus  available  for  the  generation  and  detection  of 
elastic  waves  in  this  laboratory.  The  equipment  which  has  been 
employed  in  this  work  was  designed  for  use  with  piezoelectric 
transducers,  which  have  normally  a  high  electrical  impedance. 

In  trying  to  generate  elastic  waves  magnetostrictiveJVf  however  , 
it  is  necessary  to  excite  the  surface  oi‘  the  transducer  with  an 
rf  magnetic  field.  This  type  of  excitation  requires  a  low- 
impedance  resonant  circuit,  and,  consequently,  considerable  effort 
has  been  directed  toward  the  modification  of  the  existing  equip¬ 
ment  to  permit  the  use  of  low- impedance  circuitry.  This  effort 
has  not  been  completely  successful  because  the  equipment  was 
designed  in  such  a  manner  that  its  modification  would  require 
a  complete  redesign.  Consequently,  a  rather  unsatisfactory 
compromise  modification  was  employed  which,  although  it  permitted 
observation  of  the  desired  phenomena,  did  not  permit  optimum 
execution  of  the  experiment.  As  a  result  of  this  inadequacy 
of  the  available  apparatus,  the  generation  of  elastic  waves  and 
their  subsequent  detection  was  rather  inefficient.  Thus,  signal- 
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tc-noise  ratios  were  indadequate  for  quantitative  measure¬ 
ments,  and  only  the  very  qualitative  results  described  above 
could  be  obtained.  New  equipment  has  been  ordered,  and  it  will 
bi  put  into  operation  during  the  second  year  of  this  program. 

This  new  apparatus  will  permit  the  quantitative  measurement  of 
elas1' ic-wave  amplitudes,  and  it  will  permit  reliable  determina¬ 
tion  of  the  polarization  (longitudinal  or  transverse)  of  the 
waves. 

It  is  clear  from  the  limited  data  which  have  been  ob¬ 
tained  that  single-crystal  transducers  of  the  rare  earths  can, 
as  expected,  be  used  as  magneiostrictive  elastic-wave  transducers. 
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Vo  THEORY  OF  MAGNETIC  ANb  MAGNETOELASTIC  RESONANCE  IN  THE 
RARE  EARTHS 

Many  attempts  to  explain  magnetic  and  maanetoelastio 
phenomena,  particularly  dynamical  phenomena  such  as  magnetic 
resonance,  in  terms  of  theories  which  had  been  developed  for 
and  applied  to  iron-group  magnetic  materials  have  not  met  with 
complete  success..  Even  the  calculation  of  some  of  the  equili¬ 
brium  characteristics  of  the  rare  earths,  such  as  the  tempera¬ 
ture  dependence  of  the  interlayer  turn  angle  in  the  helimangetis 
phase,  has  not  proved  successful  when  the  special  properties  of 
the  rare  earths  are  not  taken  properly  into  account .  As  an 
example,  calculations  of  the  field  and  temperature  dependences 

of  uniform-precession  ferromagnetic  resonance  in  dysprosium  by 
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Cooper  and  Elliott  are  not  at  ail  in  agreement  with  the  observed 
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resonance  characteristics  reported  by  Biackstead  and  Donoho 
Nor  do  the  predictions  of  Cooper  and  Elliott  agree  at  all  with 
the  very  sharp,  strongly  angular-dependent,  resonances  observed 
in  the  helimagnetic  phase  by  Biackstead  and  Donoho.  In  the  case 
of  elastic-wave  propagation  and  generation,  it  is  virtually 
impossible  to  use  existing  theories  to  explain  any  of  the  pheno¬ 
mena  observed  in  the  work  reported  here,  except,  as  Jberved  in 
Section  III,  for  the  model-independent  features  of  the  elastic- 
wave  velocity  and  attenuation  near  a  magnetic  phase  transition. 

In  an  attempt  to  explain  some  of  the  observations  made 
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in  the  course  of  the  experimental  work  reported  here,  it  became 
evident  a  major  approximation  which  has  been  made  in  most  calcu¬ 
lations  of  dynamic  magnetic  properties  in  the  iron-group  materi¬ 
als  is  totally  inappropriate  in  the  case  of  the  rare  earths » 

This  approximation,  in  which  the  energy  associated  with  the  mag¬ 
netocrystalline  anisotropy  is  approximated  in  terms  of  an  effec¬ 
tive  internal  magnetic  field,  is  valid  only  at  temperatures 
well  below  the  Curie  temperature,  and  it  is  valid  only  when  the 
anisotropy  energy  is  a  very  small  fraction  of  the  exchange  ener¬ 
gy  o  Consequently,  a  new  theoretical  treatment  of  the  magnetic 
and  magnetoelastic  properties  of  the  rare  earths  is  under  devel¬ 
opment,  and,  although  detailed  calculations  have  not  yet  been 
completely  carried  out,  the  preliminary  results  obtained  from 
this  treatment  indicate  that  several  phenomena  previously  unex¬ 
plained  or  poorly  explained  by  earlier  theories  will  be  in  good 
agreement  with  the  predictions  of  this  new  and  more  accurate 
theoretical  analysis  of  these  problems. 

Since  the  theory  is  still  under  development,  only  the 
broadly  applicable  qualitative  features  are  discussed  here* 
However,  from  the  limited  discussion  which  follows,  the  differ¬ 
ences  between  the  present  treatment  and  the  treatments  which 
preceded  it  will  be  apparent,  and  the  range  of  applicability 
of  the  new  theory  faiH  be  evident.  In  particular,  the  applica¬ 
bility  of  the  new  theory  to  the  following  types  of  problem, 


40. 


which  are  pertinent  to  the  experimental  part  of  this  program, 
is  emphasized  in  this  report: 

lo  Equilibrium  magnetic  and  magnetoelastic  properties: 

By  treating  all  the  contributions  to  the  magnetic  and 
elastic  free  energy  of  a  rare-earth  material  simultaneously, 
regarding  the  usual  magnetic  contributions,  the  exhange  and 
Zeeman  energies,  on  the  same  level  as  the  crystalline  anisotropy 
and  the  magnetoelastic  energy,  it  is  possible  to  calculate  the 
equilibrium  magnetization,  the  phase-transition  temperatures, 
the  equilibrium  magnetostriction,  and  the  helimagnetic  interlayer 
turn  angle  to  greater  accuracy  than  was  possible  in  earlier  work. 
The  reason  for  the  better  agreement  which  car  be  expected  between 
the  theory  and  the  experimentally  determined  equilibrium  values 
of  these  quantities  is  simply  that  the  earlier  theoretical 
treatments  regarded  such  quantities  as  the  crystalline  aniso¬ 
tropy  and  the  magnetoelastic  coupling  as  small  perturbations  to 
the  total  energy,  which  was  regarded  as  primarily  due  to  the 
exchange  and  Zeeman  energies,  yet,  in  the  rare  earths,  t'.  .e 
uniaxial  anisotropy  energy  is  as  large  as  one-half  the  value  of 
the  exchange  energy,  and  the  magnetoelastic  energy  is,  in  some 
cases,  ao  large  as  the  anisotropy  energy.  Perturbation  treatment 

of  these  quantities  can,  as  shown  below.  Is.  d  to  serious  mistakes 
in  the  calculation  of  equilibrium  properties. 
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2*  Magnetic  resonance  and  spin-wave  resonance: 

All  earlier  theoretical  treatments  of  spin-wave  disper¬ 
sion  relations  and  the  field  and  temperature  dependences  of 
resonance  conditions  have  made  what  appears  to  be  a  seriously 
invalid  approximation  in  the  case  of  the  rare  earths.  This 
approximation  consists  of  the  replacement  of  a  single-ion  type 
of  energy  such  as  the  uniaxial  anisotropy  by  a  term  which  is 
equivalent  to  the  Zeeman  energy  of  an  ion  in  an  effective 
field  chosen  to  give  an  energy  of  the  same  value *  This  approx¬ 
imation  is  probably  valid  in  the  case  of  small  anisotropy  at 
temperatures  far  below  the  Curie  temperature  in  a  ferromagnetic 
material,  where  each  ionic  magnetic  moment  is  strongly  aligned 
with  the  net  magnetization  by  the  strong  exchange  interaction* 

At  higher  temperatures,  however,  the  behavior  of  an  ion  subjected 
to  the  crystalline  electric  field  responsible  for  the  anistotrop^ 
is  totally  unlike  that  of  an  ion  placed  in  a  magnetic  field*  As 
a  result  of  this  widely  used  approximation,  the  calculated  spin- 
wave  properties  of  rare-earth  magnetic  materials,  which  depend 
strongly  on  the  nature  of  the  rather  large  anisotropy  found 
in  these  materials,  are  probably  considerably  in  error*  The  pre¬ 
vious  statement  was  not  more  strongly  made  because  few  spin-wave 
measurements  have  as  yet  been  reported;  comparisons  with  existing 
theory  arc,  therefore,  limited  at  this  time.  It  is  not,  however, 
difficult  to  see,  as  shown  in  whcit  follows,  that  a  more  exact 
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which  does  not  make  the  "effective-field"  approximation,  v/ill  lead 
to  results  that  are  quite  different  from  the  results  which  are 
usually  obtained  using  this  approximation. 

3.  Elastic-wave  generation  and  propagation: 

Most  previous  theories  have  generally  made  only  the 
crudest  approximations  in  order  to  find  the  coupling  between  spin 
waves  and  elastic  waves.  Normally,  the  propagation  cnaracteristics 
of  each  type  of  wave  are  calculated  separately,  and  a  coupling  inter¬ 
action,  usually  phenomenological  in  nature,  is  introduced,  leading 
to  the  usual  coupled-wave  phenomena  which  are  common  to  all  types 
of  wave  coupling.  In  particular,  the  effects  of  coupling  to  the 
magnetic  part  of  a  crystalline  system  on  the  propagation  of 
elastic  waves  have  usually  been  treated  in  terms  of  the  scattering 
of  these  waves  by  a  weak  potential  whose  effects  were  strong  only 
near  a  magnetic  phase  transition  (cf .  Section  III) .  Yet,  in  the 
case  of  the  rare  earths,  which  are  strongly  magnetos trictive,  there 
are  strong  coupling  mechanisms  which  are  totally  unrelated  to  those 
which  have  been  previously  treated  and  which  lead  to  strong  effects 
far  from  the  transition  temperatures  (both  above  and  below) . 
Conventional  theories  of  the  dynamical  magnetic,  elastic,  and 
raagnetoelastic  properties  of  solids  do  not,  however,  provide  a 
ready  quantitative  description  of  the  coupling  between  the  mag¬ 
netic  moments  and  elastic  waves.  They  cannot,  therefore,  provide 
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a  satisfactory  treatment  of  this  coupling  in  terms  of  the  known 
equilibrium  properties.  Rather,  the  usual  theoretical  approach  is 
to  treat  the  coupling  in  terms  of  a  phenomenological  interaction 
whose  relationship  to  the  known  static  properties  of  the  materials 
cannot  be  traced.  The  present  theory  attempts  to  get  around  this 
difficulty  by  deriving  the  coupled  equations  of  motion  from  a 
consistent  approach  which  leads  to  the  introduction  of  no  adjustable 
parameters  whose  sole  purpose  is  the  explanation  of  one  particular 
set  of  phenomena.  As  stated  above,  however,  the  the  ary  is  still 
under  development,  so  that  only  a  brief  description  is  presented 
here.  This  description  will  suffice,  however,  to  show  how  the 
present  treatment  leads  to  a  self-consistent  picture  of  all  the 
magnetic,  elastic,  and  magnetoelastic  properties,  both  static  and 
dynamic,  of  the  rare  earths. 

A.  Calculation  of  Equilibrium  Properties 

A  somewhat  simplified  model  is  presented  here,  for  the 
sake  of  clarity.  The  essential  features  are  retained,  however, 
and  the  model  can  readily  be  extended  to  a  more  realistic  situa¬ 
tion  characteristic  of  the  actual  rare  earths.  Here,  it  is  assumed 
that,  the  angular  momentum,  J,  per  ion  is  equal  to  1  (Planck's 
constant,  h,  is  assumed  to  be  unity  or  otherwise  incorporated  into 
other  parameters) .  It  is  further  assumed  that  the  system  of  in¬ 
terest  consists  of  an  assembly  of  such  ions,  coupled  by  exchange 
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interactions  of  the  isotropic  Heisenberg  type  and  that  each  ion  is 
subject  to  an  axial  crystalline  electric  field  which  leads  to 
an  easy  plane  of  magnetization  perpendicular  to  the  anisotropy 
axis.  The  magnetoelastic  interaction  and  elastic  energy  are 
ignored  in  this  example.  They  lead  to  no  qualitative  difference 
from  the  discussion  given  here.  The  major  point  to  be  made  here 
is  the  fact  that  the  anisotropy  associated  with  the  crystalline 
electric  field  cannot  be  treated  accurately  in  terms  of  an 
effective  magnetic  field  even  in  the  calculation  of  the  equilib¬ 
rium  magnetization  and  other  equilibrium  properties.  For  the 
purpose  of  this  example,  the  magnetization  will  be  calculated  on 
the  basis  of  the  molecular-field  model  of  the  exchange.  Thus,  the 
Hamiltonian  of  a  single  ion  becomes 

p 

3C  =  gpDrMJ  -  ~  (J  2-J  2)  +  (J  2-2/3)  ,  (4) 

B  z  2  x  y  z  ’ 

where  g  is  the  Lande  factor,  p  is  the  Bohr  magneton,  r  is  the 
molecular-field  constant,  M  is  the  magnetization,  and  is  a 
constant  representing  the  strength  of  the  crystalline  electric 
field.  The  z-axis  is  the  axis  of  spontaneous  magnetization,  per¬ 
pendicular  to  the  anisotropy  axis.  The  energy  levels  of  this 
Hamiltonian  can  be  found  exactly,  and  the  spontaneous  magnetiza¬ 
tion  can  then  be  found  by  means  of  a  straightforward  calculation 
of  the  expectation  value  of  the  magnetic  moment,  which  is  equal 
to  the  reduced  magnetization  in  magnitude.  From  this  calculation, 
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the  Curie  temperature  can  be  obtained,  together  with  the  complete 
temperature  dependence  of  the  magnetization.  One  interesting 
outcome  of  this  calculation,  taking  the  anisotropy  energy  ex¬ 
plicitly  into  account  in  the  calculation,  is  the  result  that  the 
free  energy  is  a  function  of  temperature  even  in  the  paramagnetic 
region  above  T  .  Furthermore,  the  Curie  temperature  depends 
strongly  on  the  value  of  the  anisotropy  constant,  increasing 
with  increasing  P^  until  the  anisotropy  energy  is  approximately 
one-half  the  exchange  energy.  For  large  values  of  P  ,  the  Curie 
temperature  decreases  in  this  J  =  1  case,  although  this  behavior 
is  not  significant  in  the  case  of  the  rare  earths,  which  all  have 
much  nigher  J-values.  The  effective-field  approximation  for 
the  anisotropy  replaces  the  anisot'.'py  terms  in  Eq.  (4)  by 
another  Zeeman  term,  with  an  effective  field  adjusted  to  give 
the  same  value  for  the  energy  at  T  =  0.  This  approximation,  which 
does  not  even  give  the  correct  single-ion  energies,  differs 
strongly  from  the  exact  results,  particularly  as  the  anisotropy 
is  increased. 

If  the  magnetoelastic  energy  and  the  elastic  energy  are 
included,  in  order  to  find  the  equilibrium  magnetostriction,  the 
problem,  for  J  =  1 ,  is  still  exactly  soluble  within  the  molecular- 
field  approximation  (not  to  be  confused  with  the  effective-field 
approximation,  the  molecular-field  approximation  is  much  better 
and  much  more  readily  justified  theoretically) .  In  this  case,  the 
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results  for  the  magnetostriction  in  the  easy  plane  are  qualitatively 
the  same  as  those  for  the  effective-field  approximation,  but  the 
results  for  the  axial  magnetostriction  are  quite  different,  in 
that  the  effective-field  approximation  leads  to  no  axial  magneto¬ 
striction,  whereas  the  exact  solution  predicts  an  axial  magneto¬ 
striction  comparable  to  the  planar  magnetostriction,  in  agreement 
with  experimental  observations  for  the  rare  earths. 

Work  is  in  progress  to  carry  out  detailed  calculations 
for  the  case  of  rare  earths  with  higher  j-values  and  with  all  the 
anisotropy  and  magnetoelastic  terms  appropriate  to  these  materials. 
In  principle,  this  exact  calculation  can  be  carried  out  for  all 
realistic  cases,  buc  the  numerical  complexity  increases  rapidly 
as  j  increases,  so  that  the  effort  at  present  is  being  directed 
toward  a  generalized  calculation  using  a  large-scale  computer. 

B.  Calculation  of  Dynamic  Magnetic  and  Magnetoelastic 
Properties 

A  widely  used  method  for  the  calculation  of  the  dynamic 

•t 

properties  of  a  quantum-mechanical  system  involves  the  use  of  the 
Heisenberg  equations  of  motion,  in  which  the  time  derivative  of 
the  expectation  value  of  a  quantum-mechanical  operator  is  related 
to  the  expectation  value  of  its  commutator  with  the  Hamiltonian 
operator  for  the  system  in  the  following  way: 


<P>  =  i<[K,F}> 


(5) 
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These  equations  of  motion,  when  evaluated  for  all  dynamical  variables 
in  the  system  of  interest,  can  often  be  used  to  explain  the  overall 
dynamical  behavior  of  a  system  as  it  deviates  from  its  equilibrium 
condition.  A  well-known  example  is  the  derivation  of  the  famous 
Bloch  Equations  for  nuclear  magnetic  resonance,  which  describe 
the  behavior  of  an  elementary  magnetic  dipole  in  a  magnetic  field. 

If  the  equations  of  motion  are  worked  out  for  a  system 
described  by  the  Hamiltonian  of  Eq.  (4) ,  however,  they  are  some¬ 
what  more  complicated  than  the  usual  Bloch  equations  because  of 
the  crystalline  anisotrop  terms  in  the  Hamiltonian.  For  the 
case  in  which  J  =  1,  it  is  necessary  to  evaluate  eight  equations 
of  motion  instead  of  the  three  which  would  be  needed  to  describe 
the  motion  of  a  free  dipole.  These  equations  of  motion  are  the 
following : 


J  =  -ujJ  4-  P  J 
x  y  2  - 


j 

y 


=  U)J 

X 


■  “V* 
=  0 

= 


=  2».J12+P?Jz 

■  “VVx 


-a'J 


4 


(6) 


The  new  operators  appearing  in  the  above  equations  of  motion  are 


as  follows: 
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,2  2 
J.  _  =  J  -  J 
12  x  y 


J.  a  J  -  J(J+l)/3 


J  =  J  J  +  J  J 
4  y  z  z  y 


(7) 


J_  =  J  J  +  J  J 
5  z  x  x  z 


Jr  =  J  J  +  J  J 
6  x  y  y  x 


The  abbreviation  uu  =  gu  D4  has  also  been  employed  in  the  equations 

B 

of  motion. 


It  should  be  observed  that  the  equations  of  motion  couple 
the  usual  operators,  J^,  J  ,  and  ,  to  the  new  second-degree 
operators  defined  above,  in  contrast  to  the  case  of  a  free  ion 
with  no  crystal-field  anisotropy.  The  effective-field  approximation 
mentioned  above  consists  in  the  approximation  of  the  term  involving 
in  the  equation  of  motion  for  by  the  term  2< J^/ J^.  With 
this  approximation,  the  first  of  Eqs.  (6)  becomes 


J  =  - (uj-2P  < J  >)J  ,  (8) 

x  2  z'  y  ’ 

and  the  equations  of  motion  for  J  and  J  are  no  longer  coupled  to 

x  y 

the  remaining  equations.  They  can  be  solved,  in  this  approximation, 
for  the  normal  frequencies  of  resonance,  just  as  in  the  case  of 
the  free  dipole.  These  frequencies  differ,  however,  from  the 
correct  frequencies  which  would  be  found  from  the  solution  of  the 
complete  set  of  equations  of  motion,  and  the  difference  between 


the  two  solutions  becomes  quite  large  when  the  anisotropy  constant, 
P2>  is  appreciable  in  comparison  with  the  Zeeman  frequency,  uj. 

The  exact  solution  of  the  equations  of  motion  given 
above  is  in  agreement  with  other  exact  methods  of  solution  (such 
as,  for  example,  the  diagonali2ation  of  the  Hamilconian) .  In 
this  respect,  the  equation -of -motion  approach  introduces  nothing 
new,  if  only  the  resonant  frequencies  of  a  single  ion  or  paramag¬ 
netic  assembly  of  ions  are  required.  When,  however,  the  resonant 
frequencies  of  a  ferromagnetic  assembly  of  ions  is  required,  so 
that  demagnetization  effects  and  collective  effects  must  be  taken 
into  account,  the  equation -of -mo tion  approach,  without  making 
the  effective-field  approximation^  will  provide  accurate  results 
which  should  agree  well  with  experimental  observations,  when 
magnetoelastic  coupling  is  added  to  the  problem,  this  equation- 
of -mo tion  approach  makes  the  solutions  of  problems  involving  the 
propagation  or  the  generation  of  elastic  waves  quite  straightforward. 
If  the  elastic  anc  magnetoelastic  energies  are  added  to  the 
Hamiltonirn  of  Eq.  (4) ,  then  a  coupled  set  of  equations  of  motion 
involving  both  the  crystal  strain  variables  and  the  angular- 
momentum  variablss  can  be  developed.  The  solution  of  this  complete 
set  of  magnetoelas tic  equations  will  lead  to  the  solution  of 
problems  in  magnetic  resonance,  elastic-wave  propagation,  and 


elastic-wave  generation. 
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Unfortunately,  the  complete  set  of  coupled  magnetoelastic 
equations  of  motion  leads  to  a  formidable  problem  in  numerical 
analysis,  and  solutions  have  not  yet  been  obtained.  It  is,  however, 
possible  to  see,  in  a  rather  qualitative  fashion,  that  this 
approach  will  lead  to  agreement  with  experimental  results  in  all 
the  areas  of  interest  to  this  work.  Results  are  expected  early 
in  the  second  year  of  this  program. 
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VI.  SUMMARY  AND  CONCLUSIONS 

As  indicated  in  the  preceding  sections  of  this  report, 
substantial  progress  has  been  made  in  several  areas  of  research 
related  to  the  magnetic  and  magnetoelastic  properties  of  rare- 
earth  materials.  The  studies  of  elastic-wave  propagation  and 
generation  have  revealed  thac  the  dynamic  magnetoelastic  inter¬ 
actions  are  strong  at  high  frequency  as  well  as  at  low  frequency. 
Single-crystal  magnetostrictive  transducers  have  been  constructed 
and  successfully  used  for  the  generation  of  high-frequency  elastic 
waves.  Progress  has  been  made  on  a  new  theory  which  will  explain 
most  of  the  phenomena  which  have  been  observed  experimentally. 
Finally,  the  capability  for  growing  large  single  crystals  of 
various  rare-earth  materials  has  been  developed.  With  this 
facility,  it  will  be  possible  to  study  the  properties  of  a  much 
wider  range  of  materials. 

It  is  expected  that  the  work  carried  out  in  this  program 
will  have  a  bearing  on  the  development  of  new  rare-earth  permanent 
magnet  materials  in  other  programs  sponsored  by  ARPA.  All  of  the 
materials  currently  of  interest  as  permanent-magnet  materials 
exhibit  strong  magnetostriction.  This  magnetostriction  will  un¬ 
doubtedly  have  an  effect  on  the  performance  of  magnets  fabricated 
from  these  materials,  and  it  is  hoped  that  the  fundamental  studies 
of  magnetostriction  being  carried  out  as  a  part  of  this  program 
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will  yield  information  of  value  in  the  design  of  permanent- 
magnet  devices  based  on  these  new  materials. 


/ 
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Fig.  3  -  Schematic  diagram  of  the  experimental  arrangement  for 
”  measuring  ultrasonic  velocities. 
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Temperature  dependence  of  the  adiabatic  elastic  constant 
in  terbium.  «  determined  from  a  single  measurement; 
©  determined  indirectly  from  four  measurements. 
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Temperature  dependence  of  the  parallel  (K  ),  perpendicu¬ 
lar  (K  )  and  volume  (Ky)  compressibilities  in  terbium.  The 
compressibilities  were  determined  using  the  relationships 
in  the  footnote  to  Table  II. 


of  ultrasonic  longitudinal  wave  propagating  in  the 
[0001 j  direction  near  the  Heel  temperature  in  terbium. 
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Fig.  11  -  Temperature  dependence  of  the  adiabatic  elastic  constant 
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in  holmium. 
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Temperature  dependence  of  the  directional  eonpressi- 


Dili ties  K ^  and 
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^  and  the  volume  compressibility  K^. 


of  holmium. 
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TEMPERATURE  [°K] 


Fig.  15  -  Temperature  dependence  of  ultrasonic  echo  height  in 
terbium  single  crystals  at  various  applied  magnetic 
fields .  Longitudinal  waves  au.  20  MHz  propagating 
along  the  c-axis  and  magnetic  fields  applied  in  the 
basal  plane  parallel  to  the  a-axis. 
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Fig,  16  -  Temperature  dependence  of  ultrasonic  echo  height  in 
terbium  single  crystals  at  various  applied  magnetic 
fields.  Longitudinal  waves  at  20  MHz  propagating 
along  the  c-axis  and  magnetic  fields  applied  in  the 
basal  plane  parallel  to  the  b-axis. 
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MAGNETIC  FIELD  [KOe] 


Fig.  17  -  Peak  temperature  of  the  ultrasonic*  echo  height  in 

terbium  single  crystals  (shown  in  Figs 27  and23)  versus 
magnetic  field  applied  in  the  a-  or  in  the  b-axis.  The  o 
intercept  temperature  of  the  two  straight  lines  is  228. K. 
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Temperature  dependence  of  ultrasonic  attenuation 
in  terbium  single  crystals;  longitudinal  waves  at 
20  MHz  propagating  along  the  c-axis,  without  and 
with  magnetic  field  applied  parallel  to  the  a-  or 
b-axis . 


Log-Log  plot  of  the  critical  ultrasonic  attenuation 
in  terbium  single  crystals  and  temperature.  L°r‘SJ- 
tudinal  waves  at  20  MHz  propagating  fLong  the  c- axis 
without  and  with  magnetic  field  app-ied  parallel  to 


the  a-  or  the  b-axis. 
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Pig.  20  -  Temperature  dependence  cf  the  ultrasonic  attenuation 
in  holmiur.i  single  crystals;  longitudinal  waves  at 
20  MHz  propagating  along  the  c-axis. 
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